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ABSTRACT 

This document presents the course content for a 
workshop that integrates the use of the computer algebra system 
Derive with topics in matrix and linear algebra* The first section is 
a guide to using Derive that provides information on how to write 
algebraic expressions, make graphs, save files, edit, define 
functions, differentiate expressions, calculate both definite and 
indefinite integrals, calculate Riemann sums, plot polar graphs, and 
print* Subsequent sections discuss the following topics and their 
development in Derive: (1) matrix algebra; (2) exercises in systems 
of linear equations; (3) exercises in linear algebra; (4) solutions 
to maximization problems in linear programming; (5) graphing 
strategies in linear programming with Derive; (6) examples and 
exercises with graphing; (7) the Simplex Method of solving linear 
programming problems with Derive; (8) examples and exercises with the 
Simplex Method; (9) linear programming and the Simplex Method in 
Derive; and (10) a Derive stack for LINPRO.MTH. An additional 
resource list with 16 references is provided* (MDH) 
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The Derive Workshop: Matrix Algebra and Linear Algebra was a 
workshop where participants gained knowledge of the Derive soft- 
ware with specific applications to matrix and linear algebra. 

The Derive software is a mathematical assistant for personal 
computers. It is available for purchase from 

SOFT WAREHOUSE 
3660 WAIALAE AVENUE 
SUITE 304 

HONOLULU, HAWAII, USA 96816-3236 

PHONE: (808) 734-5810 
FAX: (308) 735-1105 

Derive software is a menu driven computer algebra system 
that simplifies, solves and plots mathematical expressions. The 
system is easy to learn and provides the user with a tool for 
working with mathematics on an elementary or advanced level. It 
will perform symbolic operations in algebra, trigonometry, calcu- 
lus, vectors and matrices. The graphing package allows for 2 or 
3 dimension plots. 

The software also contains many utility files which are used 
for specialized topics such as Bessel functions. Orthogonal 
polynomials and Fesnel integrals. 

The user also has the option of creating utility files to 
deal with a topic requiring mathematical operations. This paper 
includes such a utility file designed to solve linear programming 
problems. The utility file is designed to encourage the userto 
master the techniques necessary to solve a linear programming 
problem while eliminating the arithmetic. 

This utility file on linear programming may be used or 
shared by anyone who has purchased the Derive software. It is 
important to note the Derive software is not a shareware and must 
be purchased. 
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STARTING AND STOPPING. To priuc graphs, call up graphics.com (type "graph- 
ics" (RET)) before the command "DERIVE" (RET) to start the program. Hercdes graph- 
ics cards would require the command "hercules" before startup. When you are finished, 
Save if desired and "Q" to quit ("Y" to abandon expressions). 

USING DERIVE. 

1. Commands available to you are listed at the bottom of the screen in the COMMAND 
menu. To select a conunand, you type the letter that is capitalized in the word. You 
caul also select the highlighted conunand by hitting the rettun, or enter key (RET). 
To move the highli^^ht to the desired command, hit the space bar. WHEN YOU 
PANIC: touching the escape key (ESC) several times will get you back to this original 
menu. WHEN YOU ARE DONE: the Quit command will exit Derive after asking 
you whether you intend to save the expressions you created (see the remark below on 
saving files). 

2. The conmxand Author allows you to create expressions. Hit A (for Author), then type 
in your expression. After you (RET), the expression is entered on a Ust, or stack, 
and has a line mmiber. If you have made a typing error, you can jtist Author again. 
Editing coir jiands are discussed later. 

3. The command Window allows you to manage yoiu: screen space. After selecting this 
command, you are given other commands. Next, Designate, Close, and SpUt are the 
subcommaoids we will use most. W Next moves you to the next window, £ts will the 
Fl button. To create another window, use W Split. You are then prompted as to 
whether you want to split Horizc^ntally or Vertically. Select this with the appropriate 
capital letter, and type in the line number at which you wish to split (in Text mode, 
the screen is 80 x 20). the window in you are currently active is highlighted. To toggle 
between windows, use Fl key or W Next. To close a window, be sure you are active 
in that window (toggle there) and use W Close. 

4. There is a parenthetical remark about the Help system on DERIVE. If you are vising 
a networked version, only one user can select the command Help at any given time. 
U you try for Help and get a message that this is not possible, simply abort (iype a) 
and continue with DERIVE. Try for the Help command later. 

ALGEBRA. To Author algebraic expressions, you need to know what sjrmbols to type. 
Here are some algebraic symbols and their meaning: 

1. — z minus z 

2. z + w z plus w 
Z. z — w z minrs w 

4. z ^ w OT zw z times w 

5. z/w z divided by w 



6. z Aw z reused to the power W 

7. z% z percent = z/100 

8. z\ z factorial 

Remember that the order of operations is parentheses, Ubrary functions, exponentiation, 
multiplication and division, and addition and subtraction, read left to right. Your results 
will be different if you Author 2 A 3 + 5/4 or 2 A (3 + 5)/4 or 2 A (3 + 5/4) etc. (see p.20 
of manual) 

GRAPHS* To graph an algebraic expression, it must be a function of one \'Briable. This 
is the independent variable (the inputs) and will be represented on the horizontal axis. The 
output will be represented on the vertical axis. An algebraic expression must be authored 
in an algebrziic window, and a graph must be plotted in a 2-D graphics window. The 
command Plot will automatically cause DERIVE to split the screen, with the option of 
the graph exhibited Beside, Under, or Overlaying (like a window shade) the algebra stack. 
Choose B, U, O as you desire, then DERIVE will ask you at which coltmm (row) number 
you would Uke the window to split. The given default is a halfway. See item 1 below to 
proceed from here. 

To split the screen yourself before Plotting, look above at how to split a window 
(USING DERIVE, item 3), then toggle (Fl) to the window you wotdd like to graph in. 
Use the command Window Designate to choose the 2-D designation of that window. You 
will have to abandon any algebraic expressions in that window, so make sure they are 
listed elsewhere. 

1. Graphic capabili\;y: when you start up Derive, you are in Text mode, which will only 
produce dot graphs. For something nicer, you need to move to Graphics mode. Select 
Options, then select Display. You then have several items to select. 

i. Now you need Graphics, Medium, Large, Extended, CGA. Do this by typing the 
capital letters of these options. Then when you cycle back to Graphics, RET. Your 
screen image will become larger. Some people prefer GHSEC, which keeps a smaller 
screen (more graph) but becomes monochrome. You may produce a printout of these 
graphs, see the section on printing. 

ii. Anytime you want to do lots of calctdus or matrix algebra, I suggest going back 
to Text mode (easier to read). Options Display Text High Large Extended RET. The 
command F5 (gray key) will toggle you between the two most recently chosen modes, 
usually text (what DERIVE starts in) and graphics mode. 

2. Plotting: H you are in the graphing window and choose the command Plot, Derive will 
plot whatever you have highlighted in the algebra window. If you are in the algebra 
window, you can use the arrow keys to move the highlight to the expression which 
you wish to plot. You must then choose the Plot command twice: the first P will take 
you to the plotting window and the second P will plot the expression. 

3. Removing excess plots: In the plotting window, the command Delete will allow you 
to remove some of your graphs. 

4. Scales, Cross hairs, and Zooming: The scale of a graph is noted on the bottom of a 
graphic screen, or can be foimd by pressing the Scale command (ESC after you have 
read the scale). The scale indicates the number of units per tick mark on the graph 
and may differ in the x- and y--coordiuates. The cross hair is a cross on the graph; its 
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coordinated are listed on the bottom of the screen. The move the cross, use the arrow 
keys; HOME takes the cross to the origin. Ctrl right or left arrow moves faster right 
or left, similarly with Pg Up and Pg Dn. MOVE moves the cross hair to coordinates 
you choose. CENTER centers the graph at the cross hair. To zoom out, you can press 
FIO and F9 to zoom in. See what happens to the scale as you do that. The Zoom 
command also works. 

SAVING YOUR FILES. It would be nice to save some of the work you have generated 
in a lab session, so you can pick up where you left off the next day. 

1. You need a formatted floppy disk, 5 1/4 " or 3 1/2 Insert the disk in drive B (or 
available drive) of the PC. After you have generated all the material you wish to save, 
select the command Transfer and then the command Save. You will be asked what 
kind of file (DERIVE), and the name of the file. Yes, it is true that DERIVE will 
only save lines from the stack, and not graphs or extra windows. Still, that's better 
than nothing!! Name the file B:name.mth, where name is what you want to call it. 
The file specification .mth identifies this as a file that DERIVE can read. 

2. To use a file again, after you enter DERIVE, use the command Transfer Loac^ and 
identify the name of the file (the "B:" identifies the file £is coming off your floppy). 
DERIVE will enter all the old Unes of yoxir stack. If you want to see some particular 
line from before, try the command Jtmip. It will jinnp to the line number that you 
specify. After working on your file again, you may of course save it again. The 
command Transfer Merge will append the named file to the end of an existing stack. 

EDITING COMMANDS. 

1. While you are Authoring an expression, you can edit. The Backspace key (so labeled 
or gray, with a large arrow pointing left) will delete the element to the left of the 
cursor. The Delete key (Del) will delete the character under the cursor. To move 
the cursor, use the Ctrl key simultaneously with the s or d key (one goes right, the 
other left). If you type over previous characters, they will be destroyed. To insert a 
character, you must be in insert mode-do this by pressing the Ins key. After inserting, 
go back to overwrite mode by pressing the Ins key again. Notice that the status line 
wi^' indicate whenever you are in insert mode. To Author an expression that includes 
previously existing material, try the F3 and F4 keys. While in Author, the F3 key 
will place at the cursor the material that is highlighted on the stack, F4 puts it in 
parentheses. Notice that insert versus overwrite mode will make a difference. To 
move the highlight in the stack, use the arrow keys. The right and left arrows move 
the highlight to different elements within an expression in the stack. These editing 
commands can be found in the manual on p. 11 and p. 23. 

2. Some students have akeady asked how to manage the stack better, getting rid of 
mistakes, etc. You can Remove lines from your stack, with the Remove command. 
You have to enter what line to begin the removal, TAB to the other field, and enter 
what Une to end the removal. RET, and those lines will vanish from the stack. In 
addition, the Jtnnp conmaand will jump the highlight to the line number you specify. 

FUNCTIONS. To define a function like f{x) = sin(x), you Author f{x) := sin(x). The 
tells DERIVE that you have assigned a functional name. The option Manage Sub- 
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stitute will sub in to any line niamber a value for the variable. Simplify will give a rational 
expression for the calculation, and approximate (Xj not A) will give a decimal approxima- 
tion for the calculation. To plot a function, highlight or Author just the expression, not 
the ^f{x) :="part. You have to Itansfer Clear to eliminate /(x), even if removed from the 
stack, and you will lose your entire stack, so new functions should have different names* 

DIFFERENTIATION. To differentiatie an expression, use the option Calculus Differ- 
entiate. You are asked all the important questions, such as what is to be differentiated, 
what is the variable, and what order derivative to find. After you RET, the notation is 

given, eg. — sin x on the stack, and Simplify will perform the Calculus. 

INTEGRATION. You can get DERIVE to calculate both indefinite and definite in- 
tegrals. Use the Calculus Integrate option, and either bring down the expression to be 
integrated from the stack, or backspace and delete to eliminate the line number, then type 
in your own expression (RET). The prompt will ask you for the variable of integration; 
usually this is x. DERIVE will convert this to dx that you see. After (RET), you need to 
select limits of integration; first the lower linut a (some real number) then (TAB) to type 
in the upper limit 6. If you are working on an indefinite integral, jxist hit (RET) without 
typing in any limits (don't type 0-that is a real limit). The line that appeals in the stack 
should look like integrals in the text. 

ex. Calculus, Integrate, backspace/delete the line number, sinx (RET), x , 0 (TAB) 

(alt) p (RET) will look on the stack like: / sinx dx 

Jo 

ex. Calculus, Integrate, backspace/delete the line number, sinx (RET), x (RET), (RET) 

will look on the stack like: J sinx dx 

You can ask DERIVE to Simplify either of these, although approximate is a bettor choice 
for a definite Integral. For an indefinite integral, DERIVE doesn't add +C to the encVyou 
have to \mderstsmd that it is there. 

SUMMATION* To calculate Riemann sums for large n, you use the commands Calculus 
Sum. Once again, you have many ways to enter the smnmand. Yo i can create it on a 
line, then use that line number (RET). You could backspace/delete the line number at the 
stunmation prompt and type in the expression (RET). Finally, you could backspace/delete 
the line munber, bring down a highlighted expression with the F3 key, then type in any 
additional pieces (RET). DERIVE will prompt you for an index (the "ticker") and suggest 
a possibility-type in whatever is correct (RET). Finally, you need the range of the index: 
type in the start and (TAB) to type in the final (RET). The stack will show the sigma 
notation S with your addend. You can Manage Substitute to input real nimibers for any 
of the variables in your summation, BUT DON'T EVER PUT ANYTHING IN FOR THE 
INDEX-KEEP IT "F!! Hence, if you are integrating from a to 6 and using n rectangles, 
you can Manage Substitute values for a,&, and n. To massage the sum algebraically, ask 
DERIVE to Simplify. To calculate a ntmQ.eric value (after you have substituted for a, 6, and 
n) ask DERIVE to approximate. 

ex. To calculate 1^ + 2^ + 3^ + . . . + for various n, do Calculus Sum, backspace/delete 
the line number, type in i A2 (RET), type in i (RET), type in 1 (TAB) n (RET). The 
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n 

stack should look like .Then Manage Substitute i for i (always!) and 2 for n and 

1=1 

Simplify (do you get 5?). Manage Substitute in the sum again-5 for n and Simplify. 
You can also Simplify the original simi-you get n(n + l)(2n + l)/6 which indicates 
that DERIVE knows the formxila for the sum of consecutive squares (do you?). 

POLAR AND PARAMETRIC PLOTS, Some special plots are polar coordinates 
and parametric curves, 

1. Parametric equations are Authored [x{i)^y{t)]. When you Plot, you are asked for the 
vsnge of the parameter min to max. This is how to plot a vertical line. 

2. Polztr plots require that you select Options State Polar to change the inputs to angles 
and radii. You plot functions of the angle; again you have to specify the range, min to 
max, of the angle. Remember to chemge to Options State Recteingulax to plot (x^y) 
coordinates again. 

PRINTING. The command Transfer Print Printer, executed from ein Algebra window, 
will dimip the stack to a printer. If the command "graphics" is typed prior to entering 
DERIVE, a screen dump can be made of graphs. System requirements include CGA to 
VGA graphics, no Hercules. Then hold down the S|hift and Print Screen keys simultane- 
ously to dimip the screen to a dot matrix printer. Screen dumps to a laser printer require 
a driver: such a free program is "HPPS". On a network, additional key sequences may be 
required. 



MATRIX ALGEBRA 

VARIABLE TYPE. You may choose the command Options Input Word, so that DE- 
RIVE will recognize xl as a variable (which woxild represent xi in problems). Coefficients 
will precede the variables. Be sxire that you go back to Options Input Character before 
you try something like FIT below. 

CREATING MATRICES. 

!• Use the commands Declare Matrix to enter a matrix. (Recall that you type the capital 
letter of the command you select,) You will be prompted for the dimension of the 
matrix. Type in the number of rows desired, then TAB to type in the colxunn field, 
RET. You are then prompted for the entries, proceeding along the rows. Hit RET 
after each desired entry is typed. When you have completed the process, the matrix 
will appear on the stack as you expect it to look, identified by its line nimiber. 

2. You can enter a matrix using the Author command, a good choice for small matrices, 
or vectors. Use the square braxdcets to form a vector, such as [1,2,1,2]. To create a 
matrix, you Author a vector of vectors. You list each row as a vector, and separate 
them by commas: [[1,2,1,2], [0,0,0,0]] is the matrix with 1 2 1 2 as the first row, and 
all zeroes in the second row. 

3. If you Author vector{x^^ ^)1)5), you will get a vector of elements a:'*, with n ranging 
from 1 to 5: [x, x^, x^, x^, x^]. You can generate large matrices this way if they have a 
pattern: try Authoring vector([x^ The difference between [ and ( is very 
important when you type this. There is a difference in the way you operate on vectors 
versus colmnn or row matrices. It is easier to make column matrices, especially if you 
want to transpose, etc. 

4. If you Author ideniH'!^mairix(S)^ that is what you will see on the stack. Use the 
command Simplify to see what you expect. 

MANIPULATING the STACK, If you need to use a matrix agedn, you can call it 
up several ways. While using the Author command, you can refer to a matrix by its line 
number: for exeimple, de<(#l). You can also bring dowr. the entire matrix: the F3 key 
will bring any highlighted expression on the stack down to the Author line. F4 will bring 
it down in parentheses. Experiment with moving the highlight tising the arrow keys. If 
you Author B:= and bring down a matrix with F3, then the variable B is assigned to the 
matrix in all subsequent work and may be referred to as such. 

MATRIX OPERATORS. Derive recognizes the following matrix operations. Remem- 
ber to use line nimxbers for matrix names, and after you have Authored such operations, 
the command Simplify will execute. 

1. addition is +, scalar multiplication is 2A or 2*A, cross product is cross{a^b). 

2. scalar product of vectors or product of matrices is the period . 

3. transpose is the accent ' (Notice which way it points!) 

4. exponentiation, including inverses is typed using the carat (shift 6). 

5. det[i) 

6. rovueduce{A^B) 

7. <race(A) 

8. eigenvalues{A) 

6 
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9. charpoly{Ay z) gives the characteristic polynomial of matrix A in terms of variable 
z. You co\ild use the command soLve to find the roots, which are the eigenvalues. 
Eigenvectors can be found via row reduction (see 6) or there is a command in a utiUty 
file, 

10c ahaiy) calculates the length of a vector V, 

11, element (v,n) finds the nth component of a vector v, and element(M, j, fc) finds the 
j, k entry of matrix M, Also, element(M j) finds the jth row of a matrix M, 

SOLVING LINEAR EQUATIONS. There are two ways to solve a system of linear 
equations: 

1, Form a vector (Author, and enclose in square brackets) of the linear equations. For 
example. Author [x — 3y + z = 5, x + y — 2z = 4, 2a: — 2y — z = 9], Then use the 
command soLve. 

2. Form the matrices A and h so that your system of equations is represented by Ax == 6, 
then TovjLTeduce{A^h). You must interpret your answer. 

BEST FIT. You can find the polynomial of best (least squares) fit for a given set of data 
points. Make stire you are in the defaidt variable mode, which is Options Input Character. 
Author, then Simplify, as in the example below, 

jH[[x, ax2+6x+c],[-1.5,0],[0.5,-2],[l,5,.1.5]] SimpUfies to the parabola a:V2-a:/2-^ 
You can plot the points and the approximation. 
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REMARKS* Linear Systems do not have to be square to be soLved, in either form (1-5 
or €)• 

!• Solve the following linear system: 

[a: + y = 9, 2x - y = 0] 

You Author the system as written above, using square brackets, and then soLve. Of 
course, in class one would plot the 2 lines (Author in y = form, then Plot, Plot) to 
visuahze solutions. 

2. See how DERIVE handles messy ntmibers: 

[Zx - 5y + 11.5z = 161.55, -3x + y + z = 4.7, -9x + Oy - 2z = -29.4] 

When you soLve, the solution is a rational solution; approximate gives a decimal approx- 
imation. 

3. Infinite solutions: 

[-x + y + 3Zz = 8, 5x + 2y - 7z = 37, 4x + 4y - 2z = 56] 

4. Infinite solutions with more than one parameter: DERIVE asks you which variables 
you wish dependent: 

[3x - y + llz - = -7, x + 2y-z + w = 0] 

5. No solutions: 

[x + y — z = 2, 2x — y + 3z = —5, — x + 2y — 5z := 2, x + y + z = 5] 

6. Solve the above systems using row reduction: begin with Declare Matrix. The matrix 
which you form is the augmented matrix [A\b] for the linear system Ax = 6. Of course, 
there is no "bar" in DERIVE; you have to imagine where it would be. For problem 1, 
you would Declare Matrix which is 2 TAB 3 RET, and list the entries as prompted, 
proceeding row by row. The result on line 13 (or so) of yoiu: stack is: 

1 1 9" 

2 ~1 0 

Next, you Author TOwjrtduce{# 13) (RET\ or whatever the line number is of your 
matrix. Then Simplify (RET) to obtain the reduced row echelon form. Interpretation 
of yotir results is required, but little or no back-substitution. 
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BEST FIT* FIT is an operator which finds the best fit (least stun of squares of dis- 
tances) of a curve to data points. If there axe n data points, an exactly fitting polynomial 
of degree n — 1 will be found, or best fit of polynonaials of lesser degree, FIT is Authored 
as listed in the exercised below, then Simplify or approximate. The data points can be 
plotted: move the highlight with the arrow keys up to the FIT line, then right arrow to 
highlight the matrix. When you Plot Plot, DERIVE will beep and say "too many vari- 
ables". Then DERIVE plots the data points and, imfortxmately, the y = x line. Use 
the command Delete First to remove the y = x line. Return to the Algebra window and 
highlight, then Plot the curve which represents the Fit," 

REMARK. MaJce sxire you are in the default variable type which is Options Input Char- 
acter (RET), so that DERIVE recognizes ax below as the product of two variables a and 
X, not a single variable string, 

1, Author fit[[x,ax + 6], [2, -2], [-1, 1,5], [6.2,3]] (RET). Then use the right arrow to 
highlight just the matrix part of the line etnd Plot,Plot, Tate out the y = x by Delete 
First, adjust the Scale to a; = 2 TAB y = 2 (RET) to see all the points, and return 
to Algebra window. Use the up arrow to highlight the entire Fit statement agsdn. 
Simplify the Fit., and Plot Plot, Compare to the fit with just the first two points. 
When you are finished with your plots, Delete All before returning to the Algebra 
window, 

2, Author fit[[x,ax^ + bx + c], [2, -2], [-1,1,5], [1,5], [0,37/6]] (RET), Then Plot the 
points. Algebra, Simplify and plot the fit. When you Plot the parabola, the Scale 
should be a: = 1 TAB y = 2. 

REMARKS. DERIVE is using the matrix algebra solution which " ^ves inverting a 
matrix, so the y— coordinates must be a function of x, 

3, Refer to the graph of sunrise/simset times on the last page this booklet. Gather some 
data to create a reasonalble approximating polynomial. 
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REMARKS. To create matrices, see page 6. The Declare Matrix commzmd is easiest, 
although you can Author. Every matrix is a column vector of row vectors, and vectors axe 
typed with square brackets, elements separated by commas. The Declare vectoR command 
is good for Calctilus 3, but not for Linear Algebra applications, ein n by 1 matrix works 
better (for example when you transpose). 

Use the following matrices and vectors below. Create them with Declare Matrix. 
Refer to them in calculations by line ntmiber (# 14, for example) or you caxi name them 
by Authoring A := and then using the gray F3 key to bring the appropriate matrix down 
to the author line. Then the matrix can be referred to by variable name. 



A.= 



1-14 
3 2-1 
2 1-1 



9 := 



3 .5 .45 
-4 1.5 -1 
0 0-1 



C:= 



2 TT V5 

1 -.5 2/3 



Use pi for tt amd sqrt 2 for 

v:=['2 -3 1] u:= 
1. Author Eind Simplify A + AB . 



-1 
4 

1/2 



2. Dot and matrix products are typed with a period. C.A and A.B . Compare approx- 
imate to Simplify. 

. Cu , u.v, v.u . 

4. Powers of square matrices formed by exponentiation: find A^^. The carat A is used 
for exponentiation, the shift 6 key. 

5. Author A.[[x],[y], [2:]] = u , Simplify and soLve. Notice that you retutn to linear 
systems! 

6. Check dei{A) to verify that A is invertible, then Author and Simplify A"^. 

7. Another way to solve problem 5 is to multiply the equation by the matrix representing 
A"^ To Author this, highlight the matrix that is A"^ , use F3 to bring it down without 
retyping, then type the dot product . and move the highlight up to the unsimplified 
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equation of problem 5. F4 will bring it down in parentheses. RET and Simplify. 
(Remind students to check determinaiit first; imique solutions exist if and only if the 
determinant is nonzero.) 

8. The transpose is the "infix'' operator, it is the left accent * key. Compare (A""-^)* 
to (A*)"^. Also transpose the vector v. Note that if you have commas, the vector 
doesn't transpose. For this reason, I always create vectors for Uneax algebra by Declare 
Matrix. 

9. The cross product is croaa. Here, since this is a calculus operation, you need vectors 
with commas! Find cross{\-2, -3, 1], [-1,4, 1/2]). 

10. Left multiply A by 2identHyjniatTix{Z) and SimpUfy. Experiment by Authoring other 
elementary matrices (the identity matrix changed by a single row operation) and 
multiplying them to perform row reduction. 

11. Length of a vector: Author and SimpUfy 

12. To see if the vector u is in the span of the cohimns of A, Author rowureduce (A,u). You 
can also create a new matrix, the augmented matrix, instead (see problem 6, page 8), 
as row^educe works with either one or two argtmients. The result is more than yes or 
no, you also get the coefficients! 

13. Similarly, you can use rovurfiduce to determine whether a set of vectors is dependent, 
to find the dimension of a set of vectors, and to find a change of basis matrix for 2 
bases such as the coliunns of A and B. 

14. Define a linear transformation T: here the vectors need conamas! Use Declare Func- 
tion, name it T, give it a Value [2x - 3y, -2j2r, x + 2y] . Then Author and Simplify 
T([1,0,0]), T([0,1,0]), and T([0,0,1]). Create the matrix Tl which has columns the 
image of the 3 canonical basis vectors. Verify that T([2, —9, 1.5]) yields the same value 
as the product of Tl and the vector [2 - 91.5]*. (Watch the commas!) 

15. The characteristic polynomial is foimd by Author and Simplify charpoly{A). Then 
soLve the polynomial to find its roots, the eigenvalues. You can also Author and 
Simplify eigenvalues(A). Find bases for eigenspaces by rovurtduce of the appropri- 
ate augmented matrix. By loading the utility file "vector.mth", you can find exacts 
eigenvecior{A^l) , a basis for the eigenspace of the matrix A for the eigenvalue 1 (or 
whatever A). Approx^igenvecior is appropriate for large nxuneric matrices, yielding 
an approximate eigenvector of length 1. 

16. The applications of Markov chains work nicely on DERIVE. Recall that a Markov 
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chain is represented by a transition matrix M whose ij entry represents the probability 
of transition from state i to state j. Regular transition matrices (some power of M 
contains all positive entries) have unique fixed vectors. It is easy to approximate large 
powers of M on DERIVE, and the rows of will converge to the fixed vector. One 
can also Declare Variable s with value [a:,y,2r], and Author s.M = s. Then Simplify, 
include the equation x + y + z = 1, and soLve, approX. For example, consider weather 
prediction. K the weather changes tomorrow to good, indifferent, or bad depending on 
whether it is good, indifferent, or bad today, this can be represented by a transition 
matrix, say 



M 



\ G I B 

G .5 .3 .2 

/ .3 .4 .3 

VB .2 A A 



Decide what the long term weather probabiUties are. (Don't use G,I,B in your matrix, 
of course.) 
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Solving Maximizatioa Problami in Linear Prograxaming 



The business majors at Illinois Benedictine College take a two 
semester course in finite mathematics^ which includes an 
introduction to linear prograimning. The topic inclx;des graphing 
linear inequalities in two variables, solving linear programming 
problems graphically and solving linear programming problems by the 
Simplex Method. 

Derive is an excellent tool for the graphical solution of systems 
of linear inequalities. The student quickly grasps the concept and 
has a good understanding of the region of feasible solutions and 
the objective function. 

When it comes to solving linear programming problems by the Simplex 
Method there are software packages that will solve the problem. 
Since our purpose in this course is to teach students how to solve 
linear programming problems, these packages are not suitable. 

This presentation includes a user created utility file for Derive 
that requires the student to master the strategy necessary to solve 
linear programming problems by the Simplex Method. The idea is to 
stress the technique and eliminate the arithmetic. 



1. Margaret A. Lial, Charles D. Miller and Thomas W. Hungerford, 
Mathematics with Applications in tha Managamant, Natural and 
Social Sciancas, Fifth edition; Harper Collins Publishers; 
C.1991 
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LZMKAR PROGRAMMING WITH DKRIVS 

STRATEGY WHEN GRAPHING 



•Enter the constraints as equations. 
•Solve each equation for y. 
•Decide on a feasible scale. 

•Prepare the graphing screen to display only the first quadrant. 
•Plot the lines. 

•Determine which intersections are needed. 

•Use Derive [Eq.l/ Eq.2] to find the intersections. 

•Assign two or three values to Z in the objective function and 
plot these lines. 

•Visually determine where the maximum is. 
•Calculate the maximum. 




ERIC 



14 




EXAMPLE Solution by Graphing 

Maximize: Z — Zx + 2y 
Subject to: ^ > 0 , y > 0 
a; + 2y < 120 
3x + 4?/ < 260 
2x + y< 115 

Solution: z = 190 wiien a; = 40 and y = 35. 





EXERCISES 



Maximize: z = 4x + 3y 

Subject to: x > 0 , y > 0 
y <25 

2x + 5y < 155 
3x + 2y < 150 

Solution z = 205 when x = 40 y = 15 



Maximize: z •= 2x + 2y 




Subject to: x > 0 , y > 0 
-X + 5y < 625 



« + y < 275 
3x - y < 625 

Solution z = 550 when 125 < x < 225 on the line x + y = 275 



Maximize: z = 2x + 3y 

Subject to: x > 0 , y > 0 
X + 5y < 30 
3x + 7y < 46 
X + y < 10 
2x + y < 18 

Solution z = 24 when x = 6 y = 4 
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IiZMKAR PR0GSAMMZN6 WITH DSRIVE 

STRATEGY WITH THE SII4PLEX METHOD 



•Prepare the utility file. 
Transfer, Load, linpro.mth 

•Enter constraints and objective function. 
Declare, Matrix 

Rows = # of constraints plus 1. 
Colurons = # of variables plus 1. 

•Create initial simplex plateau. 
Slackjmatrix (v^ ir j) 

•Select the pivot element. 

Most negative entry in indicator row. 
Smallest ratio, constant/coefficient . 
Divide^elements (v^ i) 

•Perform the pivot operation. 

Pivot (Vrirj) 

•Read the solution or repeat the process. 

•Solution found when indicator row contains only non-negative 
entries. 

•To facilitate reading solution. 
Dxvide_row(Vrir j) 
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EXAMPLE Solution by Simplex 

Maximize: Z — Zx\ + "^^i 

Subject to: a;i > 0 , 0:2 > 0 
a?: + 1x2 < 120 
3a;i + 4x2 < 260 
2x1 + a;2 < 115 

Solution: z = 190 when xi = 40, X2 = 35, 0:3 = 10, X4 = 0, 
xs = 0 




1 
3 
2 
-3 



2 
4 
1 
2 



1 0 
0 1 
0 0 
0 0 



0 0 120 

0 0 260 

1 0 115 
0 1 0 
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EXERCISES 

Maxiinize: Z = 5xi + + 2x3 

Subject to: Xi > 0 , xj > 0 X3 > 0 
xi + 3x2 + 7x3 < 15 
2xi + X2 + 5x3 < 12 

Solution: z = ~ when Xi = y> ^2 = T' X3 = 0, X4 = 0 and Xs = 0. 



Manmize: Z = 3xi + Xj + 2x3 

Subject to: Xi > 0 , xj > 0 X3 > 0 
xi + 2x2 + X3 < 80 
3xi + 6x2 + 2x3 < 20 
3xi + X2 + 3x3 < 60 

Solution: z = 20 whenxi = ^, xa = 0, X3 = 0, X4 = xs = 0 and xe = 40. 



Maximize: Z = 2xi + 3x2 + 5x3 + 3x4 

Subject to: xi > 0 , X2 > 0 X3 > 0 X4 > 0 
4xi + 2x2 + 3x3 + 2x4 < 12 
2xi + 4x2 + X3 + X4 < 16 
8x1 + 3x2 + 6x3 + X4 < 24 
xi + X2 + 2x3 + 4x4 ^ 32 

Solution: z = 20 when Xi = 0, X2 = 0, X3 = 4, X4 = 0, xj = 0, xg = 12, 
X7 = 0 and xg = 24. 
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LIMKAR PROGRMOaNG AHD THB SIMPLEX METHOD ON DERIVE 



Derive can be used to solve linear programming problems by the 
simplex method. The commands included in this utility file are 
desig;ied to reinforce the techniques necessary to obtain the 
solution manually . The utility file contains commands that will 
solve maximization or minimization problems where there are no 
mixed constraints. 

To begin^ enter the system of inequalities as a matrix without the 
slack variables (Use Declare Matrix) . To prepare for a 
maximization problem enter the last (indicator) row with its minus 
signs. To prepare for a minimization problem enter minus signs in 
the last (constant) column and use Derive' s transpose command 
before creating the basic simplex plateau. 

There are six Derive commands in this utility file. The user must 
be familiar with four of them. The other two are necessary to 
create the basic simplex tableau. 



1. Craata tha basic siaqplaz plataau. 

Use the command slack^matrix (v, r, c) , where v identifies the 
matrix representing tKe original problem, r is the number of 
rows and c the number of columns in the original matrix. 

4c******* ************************ 4c 

SLACK_MATRIX (v, r, c) :=MIX_ELEMENTS (FILL_MATRIX (v, r) , r,- c) ^ 
**************************************** 



The slack_matrix command uses two commands, fill^matrix and 
mix_elements . The fill_matrix command attaches an identity 
matrix of the appropriate size to the original matrix and 
mix_elements places the identity matrix between the coefficients 
and the constants to create the basic simplex tableau. 



**************************************** 
FILL_MATRIX (v, r) :=APPEND (v % IDENTITY_MATRIX (r) ) 

MIX_ELEMENTS(v,r,c) :=VECTOR(IF (m_<c, ELEMENT (v,m_) , 

IF (m_=c+r, ELEMENT (v,c) , ELEMENT (v,m_+l) ) ) ,m_,r+c) 
**************************************** 
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Iiocmta thm pivot 

To find the pivot position, the constant in each row must be 
divided by the entry in the column that corresponds with the 
most negative number in the indicator row (last row) . Use the 
command divide_elements (v, c) , where v identifies the matrix and 
c is the niimber of the column targeted as the pivot. The 
divide_elements command divides the entire row by the entry in 
column c when that entry is greater than zero. The quotients 
will appear in the constants column and the column targeted as 
the pivot will contain a one if it was positive or its previous 
entry otherwise. The smallest positive quotient with a one in 
the pivot column directs the user to the pivot. 



★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

DIVIDE_ELEMENTS (v, c) :=VECTOR(IF (m_^DIMENSION (v) .ELEMENT (v,m_) , 
IF (ELEMENT ( v, m_, c) <=0 , ELEMENT ( v, m_) , 
ELEMENT (v,m_) /ELEMENT (v,m_, c) ) ) , m_, DIMENSION (v) ) 



TJmm thm pivot conand 

At this point the pivot position should contain a one. To 
place zeroes above and below use pivot (v^ i, j) , where v 
identifies the matrix, i is the row and j the column of the 
pivot location. 



PIVOT (V, i, j) :=VECTOR(IF (m_=i, ELEMENT (v,m_) , 

ELEMENT (v, m_) -ELEMENT (v,m_, j) /ELEMENT (v, i, j) *ELEMENT (v, i) ) , 

m_, DIMENSION (v) ) 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 



Raad a solution 

It may be necessary to divide a row by a constant to read a 
solution. In order to do this use the command divide_row (v, i, d) 
where v identifies the matrix, i is the row to be divided and d 
is the divisor. 

DIVIDE_ROW(v, i,d) : =VECTOR( IF (m_=i, ELEMENT (v,m_) /d, 
ELEMENT (v, mj ) , DIMENSION (v) ) 
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5* Adjustnants for problMis with siixad constraints 

When mixed constraints occur/ it is necessary to insert a minus 
sign to replace the slack variable with a surplus variable. The 
command chsCv^i^j) changes the sign of the element (i^ j) in 
matrix v. The command scale_element (v^ i^ s) is a Derive command 
from the vector utility file that multiplies a vector by a 
scalar* 



********************************* 
CHS (V, i, j) :=VECTOR(IF (m_/=i, ELEMENT (v,m_) , 

ELEMENT (v,m_) . SCALE_ELEMENT (IDENTITY_MATRIX (DIMENSION (vM ) / 
j^-1) ) ,m_, DIMENSION (V) ) 

**************************************** 

**************************************** 

SCALE_ELEMENT (v, i, s) :=VECTOR(IF (m_=i, 

s*ELEMENT(v,i) , ELEMENT (v,m_) ) ,m_, DIMENSION (v) ) 

**************************************** 
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DXRIVK STACK FOR LINPRO-MTH 

1: "Linear programming with Derive on a matrix called v" 

2: "MIX_ELEMENTS and FILL_ELEMENTS prepares the slack_matrix" 

3: MIX_ELEMENTS (v,.r,c) : =VECTOR( IF (m_<c, ELEMENT (v,m_) , 
IF (m_=c+r, ELEMENT (v, c) , ELEMENT (v, m_+l) ) ) , m_, r+c) 

4: FILL_MATRIX(v,r) :=APPEND(v\ IDENTITY_MATRIX (r) ) 

5: "SLACK_MATRIX (v, r, c) puts Slack variables in place, 
r= # rows, c= # columns" 

6: SLACK_MATRIX (v, r, c) :=MIX_ELEMENTS {FILL_MATRIX (v, r) , r, c) ^ 

7 : "DIVIDE_ELEMENTS (v, c) 

divides each element in row i by the entry in (i,c)" 

8 : DIVIDE_ELEMENTS (v, c) :=VECTOR(IF (m_=DIMENSION (v) , 
ELEMENT (v, m_) , IF (ELEMENT (v, m_, c) <=0 , ELEI4ENT (v, m_) , 
ELEMENT (v, m_) /ELEMENT (v, m_, c) ) ) , m_, DIMENSION (v) ) 

9: "PIVOT (v, i, j) Puts zeros above and below entry (i,j)" 

10: PIVOT(v, i, j) : =VECTOR( IF (m_=i, ELEMENT (v, m_) , 

ELEMENT (v,m_) -ELEMENT (v,m_, j) /ELEMENT (v, i, j) *ELEMENT (v, i) ) , 
m_, DIMENSION (v) ) 

11: "DIVIDE_ROW(v, i,d) divides row i by a number, d" 

12: DIVIDE_ROW(v, i,d) : =VECTOR ( IF (m_=i , ELEMENT (v,m_) /d, 
ELEMENT (v,m_) ) , m_, DIMENSION (v) ) 

13: "CHS(v, i,j) changes the sign of the element in row i, 
column j . " 

14: SCALE_ELEMENT (v, i, S) :=VECTOR(IF (m_=i, s*ELEMENT (v, i) , 
ELEMENT (v, m_) ) , m_, DIMENSION (v) ) 

15: CHS(v,i, j) : =VECTOR( IF (m_/=i, ELEMENT (V, m_) , 

ELEMENT (v, m_) . SCALE_ELEMENT (IDENTITY_MATRIX (DIMENSION (v M ) , 
j, -1) ) ,m_, DIMENSION (v) ) 
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RESOURCE LIST 
CAS and Calculus, Linear Algebra 



[I] Calculus Laboratories with Maple: A tool, not an Oracle^ W. Batildry and J. Fiedler, 
Brooks/Cole Publishing, 1991. 

[2] Uses of Technology in the Mathematics Curriculum: A Manual of Laboratory Exer- 
cises and Resources for College Teachers^ B. Evans and J. Johnson, Oklahoma State 
University, 1990. 

[3] Learning by Discovery^ A. Solow, Editor; NSF Grant Associated Colleges of the Mid- 
west and Great Lakes College Association: Preliminary Edition. 

[4] Applications of Calculus^ P. Straffin, Editor; NSF Grant Associated Colleges of the 
Midwest and Great Lakes College Association: Preliminary Edition. 

[5] Calculus Problems for a New Century^ R, Praga, Editor, NSF Grant Associated Col- 
leges of the Midwest and Great Lakes College Association: Preliminary Edition. 

[6] Problems for Student Investigation^ M. Jackson, Editor, NSF Grant Associated Col- 
leges of the Midwest and Great Lakes College Association: Preliminary Edition. 

[7] Calculus and the DERIVE Program: Experiments with the Computer^ L. Gilligan and 
J. Marquardt, Gilmar Publishing, 1990, 

[8] Calculus Applications in Engineering and Science^ S. Goldenberg and H. Greenwald, 
D.C. Heath, 1990, 

[9] Calculus Laboratories Using DERIVE C- Leinbach, Wadsworth Publishing, 1991. 
[10] Calculus Activities for Graphic Calculators, D. Pence, PWS Kent Publishing, 1990. 

[II] Graphing Calculator Laboratory Manual for Calculus, C. Beckmsuin and T. Sund- 
strom, Addison Wesley Publishing, 1990. 

[12] Proceeding of the International Conferences on Technology in Collegiate Mathematics — 
Contributed Papers, (three proceedings), Ohio State University and Addison Wesley 
Publishing. 

[13] Computer Enhancement of the Calculus, W. Boyce and E. Luchins, Rensselaer Poly- 
technic Institute. 

[14] Bringing Calculus to Life, R. Decker and J. Williams, graphing calculator and CAS 

versions. Prentice Hall, 1992. 
[15] Discovering Calculus with DERIVE, J. Johnson and B. Evans, Wiley, 1992. 
[16] Linear Algebra Experiments using the DERIVE Program, M. Salter and L. Gilligan, 

Gihnar, 1992. 
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